Dipole-exchange spin waves in Fibonacci magnetic multilayers 
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A microscopic model is employed to calculate the spectrum of spin waves in quasiperiodic mag- 
netic multilayers in the dipole-exchange regime. Results are presented for structures in which thin 
ferromagnetic films are separated by non-magnetic spacers following a Fibonacci sequence and ex- 
tend previous magnetostatic calculations. The results show the splitting of the frequency bands 
and the mode mixing caused by the dipolar interaction between the films as a function of spacer 
thickness, as well as the fractal aspect of the spectrum induced by the non-periodic aspect of the 
structure. 



I. INTRODUCTION 



The development of ultrathin film deposition tech- 
niques has allowed the production of layered materials 
with properties that are not found in bulk media. Lay- 
ered magnetic structures, in particular, have attracted 
a great deal of attention due to their possible use in 
novel microelectronic devices 0. The transport prop- 
erties, as well as the collective excitations, of magnetic 
structures in which magnetic films arc interspersed with 
non-magnetic layers were shown to be quite distinct from 
those of single films. These new properties are a con- 
sequence of the fact that the interactions between mi- 
croscopic elements of the system can be influenced by 
the structure of the multilayer. Thus, by adjusting the 
composition and the thickness of the films, one gains the 
ability to modify the spectrum of excitations of the struc- 
ture. One way of achieving that control is by changing 
the thickness of the non-magnetic layers, since that can 
alter the strength of the dipolar field between the mag- 
netic films. Several studies have investigated the effect 
of this long range interlayer interaction on the magne- 
tostatic modes of magnetic structures 0, 0, 0, Ej- Re- 
cently, a microscopic theory was developed for the regime 
in which the dynamic aspect of the exchange interaction 
becomes relevant . In this dipole-exchange regime, the 
dipolar interaction acts to couple the volume spin wave 
(SW) modes of the magnetic layers and thus causes sig- 
nificant modifications of the SW frequency branches, in 
comparison with the results for the magnetostatic case. 

So far, the microscopic dipole-exchange SW theory has 
been applied to layered media in which the thickness of 
the spacers is assumed to be constant throughout the 
structure. On the other hand, several studies have fo- 
cused on the properties of layered structures in which the 
width of the elements is varied along its thickness. The 
interest on these structures was motivated by reports on 
the existence of a quasicrystaline phase in metallic alloys 
0, which have been shown to display long range orien- 
tational order, but are not invariant under lattice trans- 
lations. One type of artificial quasiperiodic structured 
considered is the Fibonacci multilayer, i.e. a multilayer 
with elements that have widths that vary according to a 
binary substitutional sequence The elementary ex- 



citations of these multilayers have been extensively stud- 
ied, both theoretically as well as experimentally (for a 
review, see [Trj|V 

To date, work on spin excitations in quasiperiodic mag- 
netic media has been restricted to the magnetostatic 
regime, in which the spins of the magnetic layers precess 
uniformly [UIH. These calculations, therefore, do not 
address the influence of the multilayer structure on the 
standing volume SW modes of each layer, which have am- 
plitudes that vary along the thickness of the films. The 
influence of the exchange interaction is also expected to 
become particularly relevant for structures with ultrathin 
layers (i.e., with films with thickness of a few monolay- 
ers). Moreover, the dipole-exchange wave vector regime 
in thin films has recently become experimentally acces- 
sible, due to the development of new experimental tech- 
niques that allow the observation of the SW dispersion 
across the whole Brillouin zone ^| . 

In this paper we use the microscopic model presented 
in Ref. Q to the study of magnetic excitations in Fi- 
bonacci multilayers in the dipole-exchange regime. The 
description takes into account the position dependence of 
the microscopic spins in each layer. This model has been 
applied previously to obtain the spin wave spectrum of 
ultrathin films of ferromagnets 01 an d antiferromagnets 
|l5|. The paper is structured as follows: in section 2 the 
structure of the multilayer is described and a microscopic 
Hamiltonian for the system is introduced. In section 3, 
numerical results are presented and discussed for multi- 
layers of GdCbj and EuO. In section 4 the main results 
are summarized and a brief conclusion is presented. 



II. MODEL 

Let us consider a system containing Ni ferromagnetic 
films separated by non-magnetic films (spacers) , as shown 
in Fig. 1. In the following discussion, the ferromagnets 
and spacers are referred to as films, whereas the spin lay- 
ers in each ferromagnet are referred to as atomic layers. 
The films have a simple cubic crystal structure, with lat- 
tice constant a and have ideal interfaces corresponding 
to (001) crystal planes. Each of the magnetic films is 
a single domain and contains N m atomic layers (in the 
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xy plane) with localized spins, which are assumed to be 
aligned in-plane and to interact by isotropic exchange 
with its nearest-neighbors and via dipolar coupling with 
all other spins in the structure. All magnetic films have 
the same thickness and composition and every site in the 
magnetic lattice has the same number of nearest neigh- 
bors, with the exception of the uppermost and bottom 
films, which are assumed to have vacuum interfaces. 

In order to obtain a Fibonacci multilayer, we assume 
that the thickness of the spacers can have two different 
values. Thus, the non-magnetic spacer films are identi- 
fied as either A or B, according to their thickness, which 
is taken as d for the B films, whereas the A spacers have 
thickness d + S. The thickness of the magnetic films, D, 
is kept constant. The spacers are then arranged as a bi- 
nary string, according to the Fibonacci sequence. This 
sequence can be generated by joining the strings accord- 
ing to the recursive rule (n > 1): 

S n +1 — {S n Sn-l}, (1) 

where Sq = B, Si = A and n is the generation number. 
The strings can also be obtained in an equivalent way by 
means of the inflation rule 

B —> A, A — > AB. (2) 

The total number of elements A and B in each sequence 
is equal to a number F n , which is given by the recurrence 
formula F n+ i — F n + F„_i, with F — F\ — 1. In the 
limit n — > oo the ratio F n j 'F n —i approaches a character- 
istic number = (1 + v5)/2, which is the well known 
golden mean. The multilayer is constructed in such a 
way that both the upper and bottom films are magnetic, 
which means that, if the total number of spacers is F n , 
the number of ferromagnets is F n + 1 (see Fig. 1). 

The SW dispersion calculation is based on a Hamil- 
tonian formalism in which the dipolar coupling between 
localized spins is included explicitly. The Hamiltonian 
for the system is 

H = - JjjSj ■ Sj - gji B H S- 

i.j i 

+(») 2 EE< S ^!, (3) 

where Jij is the exchange between nearest-neighbor sites 
i and j in the same film Hq is the Zeeman field, taken as 
parallel to the z— direction. The last term in the Hamil- 
tonian represents the contribution of the dipolar interac- 
tion. The D™J^ are the long-range dipolar coupling coef- 
ficient between any sites I and m in the same or different 
films. The a and f3 indices denote components x, y or z; 
g is the Lande factor and [Ib is Bohr's magneton. The 
expression for the dipolar factors is 

n ap _ [\rim\ 2 S a p-3rf m rf m ] 
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1: 


Schematic depiction of 


a Fibonacci 


magnetic multi 



layer. The the black dots represent ferromagnetic spins, and 
the white dots are non-magnetic sites. In this figure, the mag- 
netic films have thickness D = 2a, the B spacers have widths 
d = 2a and the A spacers have thickness d + S, with S — 2a. 



where the vector r; m = rj — r TO connects magnetic sites in 
the lattice. In order to calculate the SW frequencies, one 
can write down the equations of motion for the operators 
in the films. These are then transformed to a represen- 
tation involving a two-dimensional in-plane wave vector 
k = (k x ,k z ) parallel to the film surfaces. The resulting 
system of equations can be solved numerically for the 
frequency u> of the modes. The Fourier amplitudes of 
the dipolar terms can be expressed in terms of rapidly 
converging summations, as shown in Ref |l4j . For multi- 
layers, the expressions for the amplitudes are modified by 
the presence of the spacers. By increasing the thickness 
of the spacers, the dipolar coupling between the ferro- 
magnetic films can eventually become negligible, and the 
resulting SW spectrum must then reproduce the results 
for the single film. The spin Hamiltonian can be rewrit- 
ten in terms of boson creation and annihilation operators 
by means of the Holstein-Primakoff transformation |l6j 
which, for low temperatures (i.e. T « T c , where T c is 
the Curie temperature of the ferromagnets) can be writ- 
ten as Si + « \/2S at, Si~ » \/2S ,Si Z — S — a^cii, 
where and at are boson creation and annihilation op- 
erators, respectively. 

Next, the equations of motion for the creation and an- 
nihilation operators can be obtained, and be transformed 
to a representation involving a two-dimensional (2D) in- 
plane wave vector k„ = (k x , k y ), where n is an index as- 
signed to each magnetic atomic layer of the system. The 
Fourier transforms of the dipole sums in the Hamiltonian 
are similar to the terms calculated for a single ferromag- 
netic film 01 1 with additional terms including extra dis- 
tance factors, which correspond to the dipolar interaction 
between spins in different films. In the present case, the 
different values of thickness of the spacers correspond to 
a quasiperiodic modulation of the distance factors. This 
contrasts with the previous theories of Fibonacci multi- 
layers, which were based on transfer matrix formalisms. 
The present method allows us to introduce the quasiperi- 
odic aspect in the structure in a more straightforward 
way, by means of a simple modification of the dipole sums 
expressions of a single film. Specifically, the expressions 
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for the dipole sums that describe the interactions between 
spins in different atomic layers contain exponential terms 
of the type 

e- 2 \y\~i°\ (5) 

where \y\ is the distance between the layers, j a b is a func- 
tion of the in-plane wavevectors and a and b are dummy 
indices ^ij. Thus, for a pair of spins located in atomic 
layers of different films, an extra distance, corresponding 
to the thickness of the spacers located between said films, 
must be added to \y\, which would otherwise correspond 
to the distance between the spins if they were located in 
the same film. Thus, when writing the expressions for 
the dipole sums, one must keep track of the distances 
between each pair of films. If one labels the films from 1 
to F n + 1, the extra separation between two spins in two 
given films with labels t and u is 

|A(*,u)| = \t-u\(D + d) + L(t,u)5. (6) 

where L(t, u) is a function that gives the number of 
A blocks found between the films t and u, and has 
the following properties: L(t,u) — L(u,t), L(t,t) = 0, 
L(t, u) — L(t, u 1 ) = L(u, u') and (for u > t), L(t, u+ 1) = 
L(t, u) +L(t+l,u+l)- L(t + 1, u). 

Having calculated the expressions for the dipole sums, 
one can write down the equations of motion for the SW 
modes in matrix form. The SW dispersion can then be 
calculated numerically as the eigenvalues of a (2Ni N m ) x 
(2Ni N m ) matrix, the generation number. 

III. NUMERICAL RESULTS 

Solutions for SW dispersion relations were calculated 
numerically for Fibonacci multilayers with parameters 
corresponding to the ferromagnet EuO (T c w 69 K) and 
GdCl3(T c w 2.2 K). The parameters for the dipolar cou- 
pling strength (given in terms of the bulk saturation mag- 
netization) and the exchange field were, for the GdCl3, 
4ttM = 0.82 T and H ex = 0.54 T, respectively, and 
AttM = 2.4 T and H ex = 38 T for the EuO, where in 
both cases gfiBH ex = 6SJ. 

Figure 2 shows the lowest SW frequency bands as a 
function of reduced wavevector, for a EuO structure with 
35 ferromagnetic thin films (n = 8), each comprising 4 
atomic layers (D = 4a), with the B type spacers hav- 
ing thickness d = 5a, and the A spacers having d = 20a 
(5 = 15a), with zero external field. The frequencies were 
calculated in the Voigt geometry (i.e. when both the 
magnetization and the field are on the plane and k y = 
0) and given in units of GHz (using 7 = 28 GHz/T). In 
the absence of the dipolar coupling, the dispersion would 
display four branches (each consisting of 35 degenerate 
modes) , with the characteristic sinusoidal dependence on 
k. These modes arise due to the quantization of the vol- 
ume modes in each magnetic film. Since we assumed that 
the ferromagnets have a simple cubic crystal structure, 
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FIG. 2: Spin wave dispersion relation for a EuO structure 
with 35 films (34 spacers), for an external field Ho = 0.0 T. 
In this case each ferromagnet has 4 atomic layers, whereas 
the spacers have thicknesses d — 20a (A) and 5a (B), where 
a is the lattice parameter of the ferromagnets. 
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FIG. 3: Spin wave dispersion relation for a GdCl3 structure 
with 22 films (21 spacers), for an external field Hq — 0.36 T. 
In this case each ferromagnet has 5 atomic layers, whereas the 
spacers have thicknesses d = 20a (A) and d — 5a (B), where 
a is the lattice parameter of the ferromagnets. 
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and also that the strength of the exchange interaction 
was the same for all film layers, the exchange- dominated 
regime does not lead to the existence of surface SW 
modes. As in the periodic case, by introducing the dipo- 
lar interaction, the volume SW modes become coupled, 
which lifts their degeneracy and causes the appearance of 
frequency bands at small wavevectors. Furthermore, the 
dipole-dipole coupling causes the appearance of surface 
SW bands, which result from a widening of the Damon- 
Eshbach modes of a single film. However, in contrast with 
the periodic multilayer results, the non-periodic aspect of 
the present structure causes the appearance of small gaps 
and the splitting of the branches into sets of subbands. 
This behavior can be observed in Fig. 2, where the inset 
shows the structure of minigaps and subbands. Accord- 
ing to previous continuum calculations for thick films, 
the number and distribution of these subbands tend to 
form a Cantor set, with a hierarchy of bands and gaps 
[To) . This Cantor set-like structure can also be observed 
in Fig. 2. A further consequence of the quasiperiodic 
structure is the appearance of isolated modes, which are 
found between subbands. One of such isolated modes 
can be identified in the inset in Fig. 2. These isolated 
branches are usually found in quasiperiodic structures 
and have been shown previously to correspond to exci- 
tations that are neither extended nor localized, but have 
instead a critical behavior, with amplitudes that display 
a power-law aspect [l0|. 

The SW dispersion for a GdCl3 structure is shown in 
Fig. 3. In this case, the multilayer contains 35 magnetic 
films and 34 spacers. The magnetic films have thickness 
D = 5a, whereas for the spacers the values d = 20a 
(A spacers) and d = 5a (B spacers) were used, for an 
external field Ho = 0.36 T. The figure shows the surface 
SW band, along with two volume bands. Not shown 
in the figure are the two higher frequency bands, which 
display the same pattern of gaps and subbands as the 
lower ones. As in the previous figure, the results show 
the formation of frequency bands at small wavevectors, 
with the appearance of small gaps inside the bands. The 
present results, however, also show a strong mixing of the 
surface SW band (which corresponds to a broadening of 
the Damon-Eshbach modes) and the lowest lying bulk 
band. A consequence of this mixing is the presence of a 
frequency gap at small wave vectors, induced by a mode 
repulsion effect. This feature has also been reported for 
the periodic multilayer case |6j. 

Figure 4 shows a plot of the lowest SW frequencies as 
a function of the thickness parameter S for a EuO mul- 
tilayer containing 35 magnetic films, for k x = 0.037r/a, 
D = 5a and zero external field. For 6 = 0, the results 
correspond to the spectrum of a periodic structure with 
d = 5a. As the thickness of the A spacers increases, the 
results show the appearance of gaps and the formation 
of subbands. A similar behavior is found for a multi- 
layer with 35 films of GdCbj, as shown in Fig. 5, for 
k x = 0.057r/a. In this case the B spacers have thickness 
d = 5a, and the ferromagnets have D = 8a. For the pe- 
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FIG. 4: Lowest spin wave frequency branches as a function 
of the spacers thickness ratio, for a EuO structure with 35 
films (34 spacers), for an external field Ho = 0.0 T and k x = 
0.037r/a. In this case each ferromagnet has 5 atomic layers 
and the B spacers have thicknesses d = 5a, where a is the 
lattice parameter of the ferromagnets. 



riodic case (6 = 0), the results correspond to 8 frequency 
bands, which arise from the 8 SW modes of the films. As 
the thickness of the A spacers increases, the figure shows 
the splitting of the four lowest bands. The same effect 
is observed in the upper bands. For larger values of 6, 
the greater distance between magnetic films separated by 
the A spacers means that the influence of the dipolar in- 
teraction across them becomes negligible, and the results 
match the SW spectrum of bilayers (separated by the B 
spacers) and single films (i.e. films located between A 
spacers). 



IV. CONCLUSIONS 

We presented the first results for the dipole-exchangc 
spin wave spectrum of quasiperiodic magnetic multilay- 
ers in which thin ferromagnetic films are separated by 
non-magnetic spacers. The thickness of the spacers was 
assumed to follow a Fibonacci binary string created by 
a substitutional rule. Results were obtained by means of 
a microscopic spin Hamiltonian that includes the short- 
range exchange interaction and the long-range dipole- 
dipole coupling between localized spins. The dipolar in- 
teraction, by coupling the magnetic films, acts to lift the 
degeneracy of the surface and standing volume modes of 
the films, thus creating frequency bands, which are dis- 
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tributcd as in a Cantor set. The spin wave dispersions 
also show the presence of isolated frequency branches 
that have been known to be associated with critical 
modes. These results are consistent with previous con- 
tinuum calculations that obtained the long-wavelength 
modes for Fibonacci multilayers. The presence of the 
volume modes, due to the position dependence of the 
spin sites, together with the dipolar interaction, allow 
the existence of mode mixing effects, which could not be 
obtained in the continuum calculation. 



FIG. 5: Spin wave frequency branches as a function of the 
spacers thickness ratio, for a GdCl3 structure with 35 films (34 
spacers), for an external field Ho = 0.36 T and k x = 0.057r/a. 
In this case each ferromagnet has 8 atomic layers and the 
B spacers have thicknesses d — 5a, where a is the lattice 
parameter of the ferromagnets. 



The authors would like to acknowledge the financial 
support of the Brazilian agency CNPq. 



[1] Ultrathin Magnetic Structures, Vols. I and II, eds. B. 
Heinrich and J. A. C. Bland, Springer- Verlag, Berlin 
(1994). 

[2] P. Griinberg and K. Mika, Phys. Rev. B 27, (5) 2955 
(1983). 

[3] P. R. Emtage and M. R. Daniel, Phys. Rev. B 29 (1), 
212 (1984). 

[4] K. Mika and P. Griinberg, Phys. Rev. B 31 (7), 4465 
(1985). 

[5] P. Griinberg, J. Appl. Phys. 57 (1), 3673 (1985). 

[6] J. Milton Pereira Jr. and R. N. Costa Filho, Eur. Phys. 

J. B 40, 137 (2004). 
[7] D. Shechtman, I. Blech, D. Gratias, and J. W. Cahn, 

Phys. Rev. Lett. 53, 1951 (1984). 
[8] R. Merlin, K. Bajema, R. Clarke, F.-Y. Juang and P.K. 

Bhattacharya, Phys. Rev. Lett. 55, 1768 (1985). 
[9] J. Todd, R. Merlin, R. Clarke, K. M. Mohanty, and J. D. 



Axe, Phys. Rev. Lett. 57, 1157 (1986). 
[10] E. L. Albuquerque and M. G. Cottam, Polaritons in Peri- 
odic and Quasiperiodic Structures, Elsevier, Netherlands 
(2005). 

[11] J. W. Feng, G. J. Jin, A. Hu, S. S. Rang, S. S. Jiang, 
and D. Feng, Phys. Rev. B 52, 15312 (1995) 

[12] D. H. A. L. Anselmo, M. G. Cottam, and E. L. Albu- 
querque, J. Appl. Phys. 85, 5774 (1999). 

[13] R. Vollmer, M. Etzkorn, P. S. Anil Kumar, H. Ibach, and 
J. Kirschner, Phys. Rev. Lett. 91, 147201 (2003). 

[14] R. N. Costa Filho, M. G. Cottam and G. A. Farias, Phys. 
Rev. B 62, 6545 (2000). 

[15] J. M. Pereira Jr. and M. G. Cottam, Phys. Rev. B 68, 
104429 (2003). 

[16] F. Keffer, Spin Waves, in "Handbuch der Physik",vol. 
XVIII/B, ed. S. Fliige, Springer- Verlag, Berlin, 1966. 



